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ABSTRACT 

From a large database of (1) 40,000 SOFIO/MDI line-of-sight magnetograms covering the passage of 
1,300 sunspot active regions across the 30°-radius central disk of the Sun, (2) a proxy of each active 
region's free magnetic energy measured from each of the active region's central-disk-passage 
magnetograms, and (3) each active region's full-disk-passage history of production of major flares and 
fast coronal mass ejections (CMEs), we find new statistical evidence that (1) there are aspects of an 
active region's magnetic field other than the free energy that are strong determinants of the active 
region's productivity of major flares and fast CMEs in the coming few days, (2) an active region's recent 
productivity of major flares, in addition to reflecting the amount of free energy in the active region, also 
reflects these other determinants of coming productivity of major eruptions, and (3) consequently, the 
knowledge of whether an active region has recently had a major flare, used in combination with the 
active region's free-energy proxy measured from a magnetogram, can greatly alter the forecast chance 
that the active region will have a major eruption in the next few days after the time of the 
magnetogram. The active-region magnetic conditions that, in addition to the free energy, are reflected 
by recent major flaring are presumably the complexity and evolution of the field. 

Key words: Sun: coronal mass ejections (CMEs) - Sun: flares - (Sun:) solar-terrestrial relations - Sun: 
surface magnetism 


1. INTRODUCTION 

The greatest disruptions of the Earth's magnetosphere and ionosphere and the greatest solar 
energetic particle (SEP) storms in Earth's sector of the heliosphere are produced by the greatest 
explosive releases of magnetic energy on the face of the Sun. These releases are from the magnetic 
fields that have the greatest free energy, the extremely contorted fields in the chromosphere and low 
corona in certain large sunspot active regions (e.g., Canfield et al 1999). The contorted field explosively 
releases a fraction of its free energy by undergoing an eruption in which reconnection divides and 
unleashes the erupting field to simultaneously explode outward and implode inward (e.g., Moore et al 
2001). The outward exploding field becomes a large fast (super-Alfvenic) coronal mass ejection (CME), a 
huge plasmoid that explodes out into the solar wind, driving a bow shock that often produces a SEP 
storm that penetrates much of the heliosphere. In this paper, a fast CME is one that has a frontal speed 
of at least 800 km s' 1 , the frontal speed of the slowest CMEs that produce appreciable SEP events (Kahler 
2001). The release of magnetic energy in the imploding field produces a major flare. In this paper, a 
major flare is one that puts out a major burst of soft X-ray emission, a burst of GOES class M or X. 

Figure 1 shows a large fast CME blasting out into the solar wind. It was produced together with an X- 
class flare and had a speed of about 1100 km s' 1 (Yashiro et al 2004). Because the source active region 
was about 20° east of central meridian, this CME was centered on lines of the solar wind's Parker-spiral 
magnetic field that passed far eastward of Earth at 1 AU. So, Earth-Moon space was magnetically 
shielded from most particles accelerated by the CME's bow shock. Even so, this CME produced a 
moderate SEP storm in Earth-Moon space (Yashiro et al 2004). In addition, the flare's X-ray and EUV 



emission strongly disturbed the ionosphere. The ionosphere is similarly disturbed by any M or X flare 
that happens anywhere on the disk, regardless of whether the eruption also produces a CME. Most M 
or X flares are produced together with a CME, and about half of these CMEs are fast CMEs. The closer 
the source active region is to disk center the more nearly Earth-centered the swath of the CME is likely 
to be, and the greater the CME's impact on the magnetosphere is likely to be. 

Because of the disruptive and destructive effects that major flares and fast CMEs can have on global 
radio transmission, and on communication satellites and power grids caught unawares, advanced 
warning of these solar events is needed to prepare for disruptions and to safeguard these vital systems. 
This requires reliable forecasting of the chance that a solar active region will produce a fast CME and/or 
major flare in the coming day or so, forecasting based on observable indicators of the magnetic field's 
readiness to explode. 

Previous work has established that proxies of an active region's free magnetic energy can be 
measured from the line-of-sight component of the magnetic field in the photosphere, provided the 
active region is within -30° of disk center, close enough to disk center that the line-of-sight component 
adequately approximates the field's vertical component (Falconer et al 2003; Schrijver 2007). These 
free-energy proxies are measures of photospheric signatures of the free energy that is in the field above 
the photosphere. The proxy measures do not have units of energy, but each increases with the active 
region's free energy. We have found that the particular free-energy proxy used in this paper has a 
success rate of about 80% as a predictor of whether an active region will produce a CME in the coming 
few days (Falconer et al 2008). The number of major flares and the number of fast CMEs produced by 
an active region in a day are each strongly correlated with the value of the active region's free-energy 
proxy measured at the start of the day. This correlation for each of these two kinds of events forecasts 
from an active region's measured free-energy proxy the chance that the active region will produce an 
event of that kind in the coming day or so (Falconer et al 2011). 

It is observed that active regions that have had a major flare in the past day or so have more major 
flares in the coming day or so than active regions that have not recently had a major flare (Zirin 1988; 
Zirin & Marquette 1991; Wheatland 2004; Welsch 2009). It is also observed that the more flare- 
productive active regions tend (1) to have a complex multipolar arrangement of opposite-polarity 
magnetic flux rather than a single-bipole arrangement, and/or (2) to be rapidly evolving via convective 
flux transport and/or flux emergence and/or flux cancelation (Svestka 1976; Zirin 1988; Zirin & 
Marquette 1991; Welsch 2009). Evidently, in addition to depending strongly on the active region's free 
energy, an active region's flare productivity also depends significantly on the complexity of the active 
region's field configuration, on the rate of evolution of the configuration, and on the occurrence of 
emergence and cancellation of flux in the active region. We therefore expect that an active region's 
recent production of major flares reflects these additional factors along with reflecting the free energy 
in the active region's magnetic field. 

In this paper, we show that for active regions having prior major flaring and mid-range free energy, an 
active region's coming-day major-flare-production chance and fast-CME-production chance forecast 
from the active region's free-energy proxy are much greater when the prior major flaring is taken into 
account instead of being ignored. This confirms that, in addition to depending strongly on the free 
energy, an active region's production of major flares and fast CMEs depends significantly on one or 
more other aspects of the active region's evolving magnetic field. 



2. DATA 


The data used in this study are from the large database that was compiled for the forecasting tool 
presented in Falconer et al (2011). This tool forecasts from an active region's free-energy proxy 
measured from a line-of-sight magnetogram the chance that the active region will have a major eruption 
during the next day or so after the time of the magnetogram. The database is from about 1,300 sunspot 
active regions that occurred during 1996-2004, the first 8 years of sunspot cycle 23. The database 
includes each active region's full-disk-passage history of production of major flares observed by GOES 
and SOHO/EIT and fast CMEs observed by SOHO/LASCO. Each active region's passage across the 30°- 
radius (0.5 R Sun ) central disk was tracked in the full-disk line-of-sight magnetograms from SOHO/MDI, 
which have a cadence of 15/day (96 min), yielding a total of about 40,000 active-region magnetograms 
for the 1,300 active regions. Along with the event-production history of each active region, the 
database has the value of the free-energy proxy measured from each of the active region's central-disk- 
passage magnetograms. 

What we call the "gradient-weighted length of strong-field neutral line" is the free-energy proxy that 
we measured from each active-region magnetogram. This free-energy proxy is denoted by l WL sg and is 
defined by l WL sg = J I VB L0S I dl, where the integral is the line integral on all "strong-field" intervals of the 
magnetrogram's neutral lines, the intervals on which the transverse component of the potential field 
computed from the magnetogram is stronger than 150 G, and | VB LO sl is the absolute value of the line- 
of-sight field's gradient measured from the magnetogram. The measurement uncertainty in the value of 
l WL sg is determined from the uncertainties in the measured magnetogram as described in Falconer et al 
(2008). 

The quantity l WL sg is a proxy measure of an active region's free magnetic energy by virtue of (1) the 
observed strong correlation between the gradient of the line-of-sight field and the shear in the 
transverse field at neutral lines in vector magnetograms of active regions on the central disk (Falconer et 
al 2003), and (2) the observed strong correlation between the amount of magnetic shear along the 
neutral line and the net electric current arching over the neutral line in bipolar active regions (Falconer 
et al 2006). Because the net electric current in a bipolar active region directly corresponds to the free 
energy, it is certainly a proxy of the active region's free energy. Flence, because in both bipolar and 
multipolar active regions l WL sg is strongly correlated with the amount of magnetic shear along the 
strong-field neutral lines, and in bipolar active regions the magnetic shear along the neutral line is 
strongly correlated with the net electric current, l WL sg is a proxy measure of free magnetic energy for all 
sunspot active regions, bipolar and multipolar (Falconer et al 2008). 

One of the active-region magnetograms from which l WL sg was measured for the database is shown in 
Figure 2. This is one of the measured magnetograms of the active region (AR 9066) from which the fast 
CME in Figure 1 exploded. It shows the active region's magnetic field about 15 hours before the start of 
that explosion. It is from a full-disk MDI magnetogram, and covers all of the active region's line-of-sight 
flux stronger than 100 G. The colored paths trace the strong-field intervals of the neutral lines. The 
strength of the gradient ( | V B L0S I ) at each point of the neutral line in these intervals is shown by the 
shade of color, the lighter the shade the stronger the gradient. In this way, Figure 2 shows that by far 
the strongest field gradients across neutral lines in this active region were in the 8 sunspot in the center 
of the active region, and hence that most of the value of l WL sg for this active region came from this big 8 
sunspot. 



3. FORECASTING FROM FREE-ENERGY PROXY ALONE 


In this Section, we present our method of obtaining, from the database for a given type of event, a 
forecasting curve which gives for an active region the event rate expected during the next day or so as a 
function of the active region's present measured free-energy proxy l WL S g- This is the method used in 
Falconer et al (2011). Here, we employ the method to obtain two forecasting curves: one for major- 
flare productivity and the other for fast-CME productivity. 

Each forecasting curve is obtained from a histogram of the observed event rate as a function of l WL sg . 
In this work, for forecasting from l WL S g alone, the histogram we use has 20 equally-populated bins, each 
bin spanning about 2,000 consecutive ordered values of l WL sg , each value measured from a separate 
active-region magnetogram. From the database, we have for each measured magnetogram the number 
of major flares and the number of fast CMEs produced by the active region in the next 24 hours after the 
time of the magnetogram. For major flares, for example, the number of next-24 hr major flares 
summed over the magnetograms in a bin and divided by the number of magnetograms in the bin is the 
observed average next-24 hr rate of production of major flares by active regions having measured values 
of l WL sg in the span of that bin. From the 20-bin histogram, we obtain a 20-point Log-Log plot of the 
observed average next-24 hr event rate versus the average value of l WL sg in each bin. In Figure 3, the 
first plot is the one obtained for coming major-flare productivity, and the second one is for coming fast- 
CME productivity. 

In Figure 3, the diamonds show the coming-day event rate for the bins in which the average l WL sg is 
10 4 G or greater, and the asterisks show the rate for the bins of smaller l WL sg . The error bar of each 
point shows the 1 -ct uncertainty in the coming-day event rate found for that bin. We estimated this 
uncertainty by using Poisson statistics (Iman & Conover 1983) and conservatively assuming that for each 
active region only magnetograms at least a day apart (about 1 in 15 of the active region's magnetorams 
from MDI) give non-redundant samples of the active region's coming-day event rate (Falconer et al 
2011). The Poison statistics take into account both the number of non-redundant magnetograms in the 
bin and the fraction of these for which the magnetogram's active region had an event in the coming 24 
hours. Because the fraction of active regions that produce major flares and/or fast CMEs decreases with 
decreasing l WL sg and approaches zero as l WL sg decreases below 10 4 G, the error bars become 
excessively large as l WL sg decreases below 10 4 G: the uncertainty becomes increasingly larger than twice 
the measured value of the coming-day event rate. For this reason, we obtain our forecasting curves by 
fitting only the points for l WL sg above 10 4 G, the diamond points. 

In each of the two plots in Figure 3, the straight line is a least-squares linear fit to the diamond points 
in the Log-Log space of the plot. That is, the fit is for the linear dependence of Log R on Log l WL sg , 
where R is the event rate in events per day. Both the fit and the uncertainty in the fit come from both 
the Log R values of the points and their uncertainties: the greater the uncertainty in Log R the smaller 
the weight of the point in the fit. In each plot, the fit's slope and the slope's 1 -ct uncertainty are given in 
the upper left corner. These fits show that for both coming major-flare productivity and coming fast- 
CME productivity, the increase in the coming event rate R with increasing free-energy proxy l WL sg is 
approximately a constant-index power law: on the least-squares line, (Log R)/(Log l WL sg ) = a, or R oc 
( l WL sg )“, where the power-law index a is the slope of the line. 

The two fit lines in Figure 3 are our forecasting curves for forecasting an active region's coming major 
flare productivity and coming fast-CME productivity. For a presently observed active region, from the 
active region's value of l WL sg [measured from a line-of-sight magnetogram (such as from the Solar 
Dynamics Observatory (SDO)) and calibrated to give the value that would be measured from an MDI 
magnetogram], the forecasting curves give the active region's expected rate of production of major 
flares and expected rate of production fast CMEs in the coming few days. For any chosen time interval t 
spanning no more than a few days from the time of the measured magnetogram, the forecast percent 



chance P that the active region will produce a major flare (fast CME) during that interval is given by P = 
(1 - e" Rt )100%, where R is the coming event rate forecast from the active region's value of l WL S g via the 
forecasting curve for major-flare (fast-CME) production (Wheatland 2001; Moon et al 2001). 


4. CORRELATION OF PRIOR FLARING WITH FREE MAGNETIC ENERGY 

As was mentioned in Section 1, observers have long noticed that active regions that have had a major 
eruption (major flare and/or fast CME) in the past day or so are much more likely to have a major 
eruption in the next day or so than those that have had no recent major eruption (e.g., Zirin & 
Marquette 1991). In this Section, we show that an active region's rate of production of major flares in 
the past 24 hours and its rate of production of major flares in the coming 24 hours are about equally 
strongly correlated with the active region's present value of the free-energy proxy l WL sg , and likewise 
for fast-CME productivity. This brings into question whether the observed strong correlation of an 
active region's coming major-eruption productivity with its recent major-eruption productivity 
significantly reflects any dependence of coming major-eruption productivity on active-region conditions 
other than the free magnetic energy. Our answer to this question is the point of this paper, and is 
presented in Section 5. 

The two plots in Figure 4 directly correspond to the two plots in Figure 3. The only difference is that 
Figure 4 shows the correlations between the free-energy proxy l WL sg and the observed event rate for 
the past 24 hours for major flares (left plot) and fast CMEs (right plot), whereas the corresponding plots 
in Figure 3 show these correlations for the coming 24 hours. For major flares, the slope of the power- 
law fit is 1.98±0.15 for prior flaring (Figure 4) and 2.05±0.16 for coming flaring (Figure 3). For fast CMEs, 
the slope is 2.0610.30 for prior productivity (Figure 4) and 1.9010.28 for coming productivity (Figure 3). 
Thus, within their 1-a uncertainties, the slope of the fits is the same for all four plots, and is about 2. 
This result shows that an active region's prior productivity and coming productivity of major eruptions 
both depend strongly on the active region's present free magnetic energy, and that the dependence is 
about equally strong for both. This suggests that an active region's prior major flaring might not be 
significantly separate from the free-energy proxy l WL sg as an indicator of coming major-eruption 
productivity. 


5. FORECASTING FROM FREE-ENERGY PROXY AND PRIOR FLARING 

In this Section, we show that, in addition to being an indicator of the free energy in an active region's 
magnetic field, an active region's recent major-flare productivity is an indicator of other conditions that 
are significant determinants of the active region's coming productivity of major eruptions. To show this, 
we divide our set of 40,000 active-region magnetograms into two subsets: the magnetograms of active 
regions that had a major flare within 24 hours before the time of the magnetogram, and the 
magnetograms of active regions that did not have a major flare in the past 24 hours. We refer to these 
two subsets as the "prior-flaring subset" and the "no-prior-flaring" subset. By the same method as we 
obtained the plots in Figure 3 from the full set of magnetograms, from each subset we obtain the 
corresponding plots of coming-day rate of production of major flares and coming-day rate of production 
of fast CMEs as functions of the free-energy proxy l WL sg . These plots are shown in Figure 5. 

In each plot in Figure 5, the points, their straight-line fit, and its slope are shown in red for the prior- 
flaring subset and in blue for the no-prior-flaring subset. As in Figure 3, the points having l WL sg > 10 4 G 
are diamonds, the points having l WL sg < 10 4 G are asterisks, and the power-law line is fit to the diamond 
points. The black line in each plot in Figure 5 is the power-law line from the corresponding plot of Figure 



3 for the full set. Because only a small fraction of all active regions produce major flares, only about 4% 
of the full set of 40,000 magnetograms is in the prior-flaring subset and about 96% are in the no-prior- 
flaring subset. So, for the plots for the no-prior-flaring subset we divided this subset into 20 equally- 
populated bins as in Figure 3 for the full set. But, because the prior-flaring subset is 25 times smaller 
than the full set, to keep the error bars for the prior-flaring subset less than about twice the event rate 
for the points having l WL sg > 10 4 G, we divided the prior-flaring subset into 10 equally-populated bins 
instead of 20. 

Not surprisingly, the Figure 5 plots for the no-prior-flaring subset show for both major-flare and fast- 
CME coming-day productivity the increase in coming-day event rate with increasing free-energy proxy 
l WL sg follows a power-law that is nearly the same as for the full set. In contrast, for each of the two 
kinds of events, the plot for the prior-flaring subset follows a power-law trend that is clearly different 
from the power-law for no-prior-flaring active regions. Over most of the range of l WL sg , the coming-day 
event rates for the points from the prior-flaring subset are significantly higher than for points of the 
same l WL sg value from the no-prior-flaring subset. Also, the slope of the power-law fit for the prior- 
flaring subset significantly flatter than that for the no-prior-flaring subset. Hence, the forecasting curve 
for prior-flaring active regions (the red line) intersects the forecasting curve for no-prior-flaring active 
regions (the blue line) near the high end of the l WL sg range, at about 10 5 G, and lies increasingly farther 
above the no-prior-flaring forecasting curve with decreasing l WL sg below 10 5 G. 

Thus, for both major-flare and fast-CME coming productivity, Figure 5 shows that active regions 
having values of the free-energy proxy l WL sg in the range from ~10 4 G to ~ 5 x 10 4 G have a much 
greater chance of producing an event in the coming day or so if they have produced a major flare in the 
past day than if they have not. This means that an active region's prior major flaring is an indicator of 
conditions in an active region other than the free magnetic energy that are significant determinants of 
the active region's coming productivity of major flares and fast CMEs. 


6. DISSCUSION 

Figure 5 shows that an active region having free energy less than the amount corresponding to a value 
of ~10 4 G for its free-energy proxy l WL sg has no more than a few percent chance of producing a major 
flare and/or fast CME in the coming day. This reflects the statistic that hardly any active regions having 
l WL sg < 10 4 G produce a major flare and/or fast CME: of the hundreds of active regions in our database 
that had l WL sg < 10 4 G, less than 1% produced a major flare and/or fast CME. We interpret this to mean 
that active regions having l WL sg ~ 10 4 G are hardly able to have a major eruption simply because they 
have hardly enough free energy, and that active regions having l WL sg much less than 10 4 G do not have 
enough free energy to have a major eruption. 

Figure 5 also shows that at l WL sg ~ 10 4 G, an active region from the small minority that have had a 
major flare in the past day are 5-10 times more likely to have a major eruption in the coming day than 
an active region from the large majority that have not had a major flare in the past day. As l WL sg 
increases above 10 4 G the difference between the two forecasting curves decreases until the two lines 
intersect just beyond l WL sg = 10 5 G. As l WL sg increases beyond its value at the intersection (at about 1.5 
x 10 5 G), the steeper expected-event-rate line for no-prior-flaring active regions rises above the flatter 
expected-event-rate line for prior-flaring active regions. In our database, for active regions having l WL sg 
above a given lower bound, the fraction of active regions that are prior-flaring active regions steadily 
increases from about 4% when the lower bound is near the bottom of the range of l WL sg (at ~ 10 3 G) to 
about 96% when the lower bound is at the intersection of the two forecasting curves. We interpret 
these results as follows. Active regions having l WL sg ~ 10 4 G or greater have enough free energy for a 



major eruption and become increasingly more easily destabilized to have a major eruption with 
increasing l WL sg and corresponding increasing free energy. Active region having l WL sg ~ 10 5 G have so 
much free energy that they can hardly keep from exploding. That is, for an active region having that 
much free energy it is unlikely that a major eruption will not be triggered in the span of a day. 
Consequently, the rare active region that has l WL sg > 1.5 x 10 5 G and has not had a major eruption in the 
past day is somewhat more likely to have a major eruption in the coming day than an active region of 
the same large l WL sg that has had a major eruption in the past day, simply because an active region's 
chance of having another major eruption increases with increasing time since its last major eruption. 

We expect that the difference in coming major-eruption rate signaled by whether an active region 
has had recent major flaring reflects the combined ability of the complexity and evolution of the active 
region's magnetic field to induce the field to erupt. Whether the complexity of the field configuration or 
the rate of evolution of the field is the more important aspect for unleashing major eruptions remains 
obscure and needs further observational study. 

This paper came from research funded by NASA's Heliophysics Division, NSF's Division of Atmospheric 
Sciences, and AFOSR's Multi-University Research Initiative. 
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Figure Captions: 


Figure 1. The large fast CME of 6 June 2000. This CME was produced in tandem with a GOES X-class 
flare by the explosion of magnetic field rooted in and around a large 8 sunspot (in NOAA Active Region 
9066) located at N20°, E18° on the central face of the Sun. Shown here is a running-difference image 
from the SOFIO/LASCO/C2 coronagraph. 

Figure 2. The strong-field intervals of the neutral lines in an MDI magnetogram of the 8-sunspot active 
region from which the CME in Figure 1 exploded. The color gives the strength of the gradient of the line- 
of-sight field ( | VB LO s I ) at these neutral lines in units of G/Mm. The lightest color is for I V B L0S I >500 
G/Mm, the range of the extreme gradients at the neutral line through the large 8 sunspot. 

Figure 3. Empirical forecasting curves for forecasting for the next day or so an active region's rate of 
production of major flares (left plot) and fast CMEs (right plot) from the active region's free-energy 
proxy l WL S g measured from a magnetogram. See the text. 

Figure 4. The correlation of prior rate of major-flare production (left plot) and prior rate of fast-CME 
production (right plot) with the active region's present value of the free-energy proxy l WL sg . The format 
is the same as in Figure 3. See the text. 

Figure 5. Curves for forecasting an active region's rate of production of major flares (left plot) and fast 
CMEs (right plot) from the active region's measured free-energy proxy l WL sg and knowledge of whether 
the active region has had a major flare in the past 24 hours. The red (blue) lines are the forecasting 
curves for prior-flaring (no-prior-flaring) active regions. See the text. 
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Main Finding 

• New evidence that, in addition to depending strongly 
on the free magnetic energy, an active region's chance 
of having a major eruption depends strongly on other 
aspects of the evolving magnetic field. 

Overview 

1 . Describe my free-energy proxy 

2. Show free-energy forecast curves for major flares 
and fast CMEs. 

3 . Show forecast curves for active regions that have 
recently produced X and M class flares. 



MSFC Vector Magnetogram of 8-Sunspot Source Region 
of a Major CME/Flare Eruption (2000 June 6) 



An active-region field’s horizontal shear is concentrated along neutral lines where the field’ 

horizontal component is strong and the vertical component’s horizontal gradient is steep. 

Observed-field upward (downward) vert. comp, is shown by solid contours or light shading (dashed 
contours or dark shading); red arrows show observed hor. comp. ; green arrows show hor. comp, of pot. 
field computed from obs. vert. comp. ; strong-observed-field (>150G) intervals of neutral lines are blue. 



Free-energy proxy from vertical-field component of vector 
magnetogram or from line-of-sight magnetogram: 

• Active-region field’s free energy is concentrated in horizontal 
shear along neutral lines invervals on which the horizontal 
compnent is strong and the vertical components horizontal 
gradients is steep. 
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Only active regions that have a large free energy are likely to 
produce major events in the next 24 hours. 
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Free Energy Only 

Recently Flaring 

Recently Non-flaring — 


Forecast of March 7 X-Flare/SEP Eruption 


• Top panel: Free Energy proxy level 
and evolution. 


• Lower Panel: Forecast M&X flare 
rate 

— Forecast using free-energy proxy only 

— Forecast using Free-energy proxy and 
previous Flare History 


• Symbols 

+ M-class Flares 
* X-class Flares & SEP 
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In Addition to Free-Energy, Something else that Persists 
on Periods of Days must be an Important Factor 


• Active Regions that are 
nonpotential today, are likely to be 
nonpotential tomorrow, as well as 
the previous day. 

• This conceivably could have been 
the sole reason that recent flare 
activity is a predictor of future flare 
activity. 



March 2012 


• That when controlling for ffee-energy proxy prior flaring still 
has additional predictive ability, indicates that some other factors 
must also play a role, a factor that also persists for days. 


Next Step: Find the additional Feactors 

Study active regions that have similar values of the 
free-energy proxy. 

Indentify which active regions are more flare 
productive and which active regions are less 
productive. 

Determined what factors, differ between the two 
subsets 

— (not total magnetic flux/active region size) 

— Evolution 

— Complexity 



Backup Slides 



Flux Content is Not an Important 
Additional Determinant 


Gray scale plot shows free 
energy/magnetic size distribution 
of 40,000 magnetograms of 1,300 
active regions. Red contours are 
0.001, 0.01, and 0.1, and 0.5 
event/day levels. 
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SRAG MAG4 Forecast Tool 

Example Display (March 6, 2012) 

• Active region in upper-left comer produced the March 7 Solar 
Energetic particle event and geo-effective CME 
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SRAG MAG4 Tool Forecast Before X5 Flare 


2012/03/06 22:23 

# AR# WL!DSG!N Lng Lat 24 Hour Event Rate Dist 

(kG) (deg) M&X CME FCME X SPE (deg) 

3 11428 8 -21 -17 0.010 0.020 0.007 0.002 0.003 27 

5 11429 69 -41 17 0.800 0.400 0.200 0.100 0.090 44! 

6 11430 14 -25 20 0.040 0.040 0.020 0.006 0.007 32! 

Disk Forecast Rates 0.900 0.500 0.200 0.100 0.100 

Multiplicative Uncertainties 2.7x 2. lx 2.3x 3. Ox 2.5x 
Disk All-Clear Forecast Probabilities 40.00% 60.00% 80.00% 90.00% 91.00% 
Uncertainties 40.00% 20.00% 10.00% 10.00% 8.00% 



